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By identifying the recently introduced Barbero-Immirzi field with the QCD axion, 
the strong CP problem can be solved through the Peccei-Quinn mechanism. A spe- 
cific energy scale for the Peccei-Quinn symmetry breaking is naturally predicted by 
this model. This provides a complete dynamical setting to evaluate the contribution 
of such an axion to the cold dark matter content of the Universe. Furthermore, a 
tight upper bound on the tensor-to-scalar ratio production of primordial gravita- 
tional waves can be fixed, representing a strong experimental test for this model. 
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According to the standard model (SM), two terms contribute to the strong CP viola- 
tion. Specifically, the non-Hermitian quark mass matrix M introduces a CP violation term 
proportional to Arg(detM). This sums up to the vacuum angle of QCD, 9, generating a 
CP violating interaction that depends on the parameter 9 = + Arg(detM). By measur- 
ing the electric dipole moment of the neutron, an extremely small upper limit can be fixed 
for 9, which turns out to be smaller than 10~ 10 . Such a small value implies an extremely 
precise compensation between two completely uncorrelated parameters: one associated with 
the global structure of the SU(3) gauge group and the other related to the SU{2) x U(l) 
breaking symmetry. The unnaturalness of this "fine tuning" goes under the name of strong 
CP problem. 

Peccei and Quinn proposed a dynamical mechanism to solve the strong CP problem pQ . 
They postulated the existence in the standard model of an additional U(l) axial symmetry, 
often denoted as U(1)pq. On the one hand, if this symmetry were exact, the CP violating 
interaction could be eliminated through a chiral rotation. On the other hand, we expect 
that the U(1)pq is spontaneously broken by the chiral anomaly. Interestingly enough, the 
Peccei-Quinn (PQ) mechanism allows us to solve the strong CP problem, even though the 
U(1)pq additional symmetry is not preserved by quantization. 

In order to briefly describe how the PQ mechanism works, it is worth recalling that 
the spontaneous breaking of the U(1)pq symmetry generates a (pseudo) Nambu-Goldstone 
boson, called axion j2], a possible cold dark matter (CDM) candidate [3J. The axion interacts 
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with matter through the following effective action [I] 
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where f a denotes the scale of the U(1)pq symmetry breaking. G = dA + ig s A A A is the 
curvature 2-form associated with the SU(3) valued connection 1-form A = A 7 A 7 , A A being 
the generators of the group, and g s the strong coupling constant. With the collective symbols 
ip and ip we denoted fermion matter fields, interacting with the axion through a derivative 
coupling term, S matt . 

The CP violating 9 term combines with the anomaly-induced interaction between the 
axion and the gluon fields; the possible observables of the theory now depend on the effective 
vacuum angle 9{x) = 9+ ?y+, sometimes referred to as the misalignment angle. The effective 
interaction, 9(x)trG A G, represents a nontrivial potential for the axion field, which selects 
a particular vacuum expectation value. In particular, the periodicity of the potential in 
the effective vacuum angle, 9(x), implies that it has a nontrivial minimum corresponding 
to 9(x) = |5J, so that (a(x)) = —f a 9. Consequently, the gluons effectively interact only 
with the physical axion aph ys (^) = o,{x) — (a(x)), preserving the theory from the strong CP 
violation [6]. 

The physical features of the axion, as, e.g., its mass and the strength of its interactions 
with ordinary matter, strictly depend on the scale of the PQ symmetry breaking, f a , which 
remains a completely free parameter, not fixed by the theory. The scope of this paper is to 
present a new model that allows us to solve the strong CP problem d la Peccei-Quinn, with 
the remarkable advantage that the parameter f a turns out to be fixed by the theory. This 
provides a completely determined dynamics, so that the contribution of such an axion field 
to CDM can be estimated as a function of the initial misalignment angle. 

Even more interesting, the model predicts the production of isocurvature fluctuations 
during inflation, allowing us to fix a tight upper limit to the tensor-to-scalar ratio, r. This 
represents an experimentally testable prediction that can, eventually, rule out the model. 

Let us start by considering a generic space-time with torsion. The chiral rotation of the 
fermionic measure in the Euclidean path-integral generates, besides the usual Pontryagin 
class, a Nieh-Yan (NY) term [7], which diverges as the square of the regulator [S] (see also 



Above, M denotes the regulator, while a is the parameter of the transformation [10]. In 
order to avoid the appearance of this divergence, one of us has recently proposed to introduce 
a field, (3(x), interacting with gravity through the Nieh-Yan density namely 
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where \ is a generic coupling constant with the dimension of energy. According to this 
proposal, we assume that ^ is the fundamental action for gravity and matter [12]. 

In order to clarify some aspects related to the proposed modification, we write below the 
resulting semiclassical effective action [TTJ Q2] : 

Sett = Sup [e] + S[A] + S mr [e, ip, W\ + \ 

+ wj * J(A) A J(A) * J{A) A ^ (4) 
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where we have defined the new field f3(x) = VGkx/3(x) and introduced the constant fz = 
H= ~ 1.98 x 10 18 GeV. An SU(3) valued connection 1-form A, representing the strong 
interaction, has been considered as well, with G being its curvature. In the last line, a trace 
over internal indexes is understood. 

Some comments are now in order. The pure gravitational sector of the effective theory is 
reminiscent of the so-called Chern-Simons modified gravity, vastly studied in the literature 
(see the interesting and complete review [H]). So, from an effective point of view, the 
modification of the gravitational action proposed in Eq. (|3| reduces to a well-known theory 
of gravity, originating from string theory and featuring some interesting dynamical effects, 
well within the presently available experimental limits [15], thus persuading us to take it 
seriously. The resulting semiclassical effective theory does not depend on the free coupling 
constant x an d shares many common features with that postulated by Peccei and Quinn. It 
is worth noting that in fact, in the case of massless fermions, the full action with the Nieh- 
Yan modification presents an additional U(1)a symmetry, which is broken at an effective 
level by the interaction between the (3(x) field and the fields strength in the last line of Q 
[16], exactly analogous to the U(1)pq introduced above. The presence of these interaction 
terms reflects the existence of the chiral anomaly, which, in fact, spontaneously breaks the 
U(1)a symmetry at the energy scale /s, naturally determined by the theory, in striking 
contrast with the Peccei- Quinn scenario, where f a is a free parameter of the theory. 

Let us now assume that the system evolves in a symmetric space-time, characterized 
by a vanishing R ab A R a b term, as the unperturbed Friedmann-Robertson- Walker (FRW) 
cosmological model. In this hypothesis, comparing action (H) with 0, one can appreciate 
the functional analogy of the two effective theories. This analogy strongly suggests to identify 
the field /3 with the axion field a (see also [17] for a supersymmetric analogous identification); 
this is the essence of our proposal, which represents the main novelty of this model, whereas, 
in some previous papers [H1EI3], the possible coexistence of the f3(x) field and the standard 
axion, a(x) was addressed. 

The term in the last line of Q represents a nontrivial potential for the Barbero-Immirzi 
(BI) axion field, selecting a particular CP preserving vacuum state. So, by implementing a 
mechanism analogous to the Peccei-Quinn one, we can solve the strong CP problem via the 
Bl-axion field. 

As was noted above, in this model the symmetry breaking energy scale fz is fixed by the 
theory, allowing us to estimate the expected zero-temperature mass of the Bl-axion field, 
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which, as for the standard axion, is generated by instantonic effects [5]. We obtain, 

/t y/ m u m d „ „, in-12 y /r-\ 

m /3 = ~r m Tr — 3.04 x 10 ev, (5) 



where we used the value of the pion decay constant, /„- = 93 MeV, measured in the decay 
process n + — >• fi + + z/ M . As is well known, instantonic effects depend on the temperature, in 
particular, we expect that the greater the temperature, the smaller the mass of the Bl-axion 
field becomes; according to the standard literature, we have that [18] 

^(£) 4 T>A, 
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H T < A, 

where we have assumed b = 0.018 and a color anomaly index equal to 1. A ~ 200 MeV is 
the QCD scale. 

So, in this model, the physical parameters, namely the mass of the field and the magnitude 
of its interaction with matter, are fixed by the theory. Remarkably, this allows us to reduce 
the parameter space of the theory and extract strong predictions from the cosmological 
scenario we are going to study. 

In general, when dealing with axion scenarios, there are two possibilities. The first one is 
that the PQ symmetry is restored after inflation, and then broken again after the Universe 
cools down. This happens if the reheating temperature Trh is larger than the energy scale 
at which the symmetry is broken. The second possibility is that the PQ symmetry is broken 
during inflation and never restored afterward. In order for the symmetry to be broken 
during inflation, the scale fpQ has to be larger than the Gibbson-Hawking temperature 
Tqh — Hj/2tt associated with the cosmological horizon (here Hj is the value of the Hubble 
expansion rate during inflation) [Tj5] ; furthermore, in order for the symmetry to stay broken 
after inflation, the reheating temperature has to be smaller than fpq. The inflationary 
expansion rate is constrained by the WMAP observations [2D] to be F 7 < 6.29 x 10 14 GeV. 
The reheating temperature is poorly constrained by the observations and could be anywhere 
in the range 1 MeV — 10 16 GeV. However in the context of the model presented here, the 
relevant energy scale fz ~ 10 18 GeV is so large that we will always have to deal with the 
second scenario, i.e., the symmetry remains broken after inflation. Nevertheless, we will also 
briefly take into account the possibility that inflation never occurred. 

The fact that the PQ symmetry stays broken after the end of inflation, has two important 
consequences. The first is that the initial misalignment angle of the BI field is practically 
constant within the region corresponding to our present horizon, and can take any value 
between — tt and n. The second is that isocurvature perturbations are produced in the BI 
field. 

The cosmological limits on axion properties have been recently reassessed in light of the 
5-year WMAP data [21] . Here we will do the same for the BI axion. Axions can be produced 
in the early Universe through two distinct mechanisms [22] , namely coherent production due 
to the initial misalignment of the axion field, and the decay of axionic strings. The latter 
is relevant only if the symmetry breaking happens after the end of inflation, then, here, we 
will be concerned only with the misalignment production. The basic idea is that, when the 
axion field is created, the initial value 9^ of the misalignment angle 9 is displaced from zero, 
since no preferred value of 9 exists. 
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Since the axion field is created during inflation, our present Hubble volume corresponds 
to a small patch at the time of creation, where the value of 9^ can be assumed to be spatially 
constant. 

In a flat FRW Universe, the zero mode of the dynamical field 9{x) evolves according to: 

i + IHi+^-O. (7) 

where a dot denotes the derivative with respect to cosmological time, H is the Hubble 
parameter, and the potential V{9) = ra~(T)/J(l — cos#) [23]. It is clear from Eq. (6) that 

in the high temperature limit T ^> A the Bl-axion is effectively massless. Then V{9) = 
and 9 = const is a solution of the equation of motion and the misalignment field is frozen 
to its initial value, 9i, until the mass becomes comparable to the expansion rate of the 
Universe, i.e. H ~ T, and the field starts oscillating around 9 = 0. For the value of the 
mass considered here, this happens at T ~ 52MeV. We have numerically integrated the 
Klein-Gordon Eq. ([7]) down to a temperature well below the onset of oscillations and used 
entropy conservation to obtain the present number density. In the limit of small 9i, this 
procedure yields 

^(To) ~ 2.8 x io^ 2 ^, (8) 
p cur 

corresponding to an energy density Pp{T ) ~ 85 9\ GeV/cm 3 . What is remarkable about this 
result is that, since the energy scale at which the symmetry breaking occurs is fixed by the 
theory, the present day energy density of the BI axion depends only on the initial misalign- 
ment angle. Given that the present critical density of the Universe is p c ~ 10~ 5 GeV/cm 3 , 
the above formula points to the necessity of having 9{ -C 1. In particular, we know from 
the recent measurements of the WMAP satellite [20J that the present dark matter density 
is Hdm/i 2 = 0.1131 ± 0.0034 at 68% CL, where Qdm = PduJ Pc is the density in units of 
the critical density, and h is the Hubble parameter in units of 100 km sec -1 Mpc -1 . Then 
assuming that Bl-axions make up all the dark matter (Qs = f2dm) the initial misalignment 
angle has to be very small: 9i ~ 1.2 x 10~ 4 . Larger values of the initial misalignment angle 
lead to a present axion density too large with respect to the WMAP value, so, in general, 
one should require that ^ < 1.2 x 10~ 4 . The axion density would be diluted, and then the 
limits on 9i relaxed, in the presence of a significant entropy production at a temperature 
below the QCD scale. This could be the case if the reheating temperature T RH < A pM] . 

Let us also examine the possibility that inflation never occurred [25]. In this case, the 
initial misalignment angle would be a function of the spatial coordinates and should be 
replaced with its average over [— n,ir], i.e. (9f) = n 2 /3. Consequently, the present axion 
energy density would be completely fixed, leading to a density parameter Q^h 2 ~ 10 8 , clearly 
overshooting the observed value by 9 orders of magnitude [26] . 

Another prediction of the model examined here is the production of axion isocurvature 
perturbations. As it happens for the inflaton field, de Sitter-induced quantum fluctuations 
in the Bl-axion field are generated during inflation. The corresponding energy density fluc- 
tuations have an amplitude proportional to Hj/(9if^) and are completely uncorrelated to 
those in the other components (radiation and matter) since axions were not in thermal 
equilibrium with photons during inflation. Moreover, since axions made a negligible contri- 
bution to the energy budget of the Universe at that time, the fluctuations in their energy 
density did not produce a corresponding perturbation in the curvature, hence the name 
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Figure 1: Constraints for the Bl-axion in the [Hi, 9i) plane. The horizontal line corresponds to 
the Bl-axion making all the dark matter in the Universe (namely 0, ~ 1.2 x 10~ 4 ). The diagonal 
line comes from the constraints on the isocurvature fluctuations [Hi/di < 8.25 x 10 13 GeV). The 
shaded area shows the allowed parameter region. The dashed and dotted diagonal lines corresponds 
to the expected improvement of the bound on isocurvature fluctuations from the measurements of 
the Planck satellite and from an ideal, cosmic variance limited CMB experiment, respectively (see 
[H] for details). 



"isocurvature" . The amplitude of primordial isocurvature perturbations can be constrained 
by observations of the CMB anisotropies; in fact, an analysis of the WMAP data yields the 
constraint H I /9 i < 4.3 x 10~ 5 //3 = 8.25 x 10 13 GeV [2T]. This bound can be combined with 
the constraint Bi < 1.2 x 10~ 4 to obtain the allowed region in the [Hi, 9j) parameter space as 
shown in Fig. [TJ In particular, the two constraints together imply Hi < 10 10 GeV. This low 
value leads to an interesting prediction of the model. Since the amount of gravitational waves 
(corresponding to tensor perturbation modes) produced during inflation is also proportional 
to Hi, the model yields an upper bound to the amplitude of tensor modes. In particular, 
in terms of the tensor-to-scalar ratio r, we get the very tight upper bound r < 1.4 x 10~ 9 . 
This means that a detection of even a very small amount of primordial gravitational waves 
by one of the upcoming CMB experiments would rule out the model proposed here. This 
would be the case, in particular, if tensor modes are detected by Planck, since it is expected 
to be sensitive to r > 0.05. 

Finally, let us briefly discuss the astrophysical constraints on the BI axion. In the case 
of the standard axion, the astrophysical constraints mainly depend on the strength of its 
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interaction with ordinary matter, since this controls the rate at which the nuclear energy 
generated in the core of a star is carried away in the form of axions, thus modifying the 
standard stellar evolution. In order to evade these constraints, the axion couplings have 
to be either small enough so that few axions are produced as a by-product of the nuclear 
reactions, or large enough in order to keep the mean free path of axions well inside the 
radius of the star. However, the computation of the couplings of the axion is completely 
general, so that the standard results also hold for the BI axion. In particular, one has that 
the couplings are inversely proportional to fz, so that we can expect them to be extremely 
small [2?]. In fact, this implies that the Bl-axion is, as long as astrophysical limits are 
concerned, roughly equivalent to a DFSZ axion with a mass m a ~ 3 x 10 -12 eV. The most 
stringent astrophysical upper limit on the axion mass comes from the observations of SN 
1987A and states that m a < 10~ 3 eV (or, equivalently, /„>6x 10 9 GeV). Thus we conclude 
that astrophysical observations cannot rule out the existence of a Bl-axion. 
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